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We present a systematic study of the microwave-induced oscillations in the magnetoresistance of 
a 2D electron gas for mixed disorder including both short-range and long-range components. The 
obtained photoconductivity tensor contains contributions of four distinct transport mechanisms. We 
show that the photoresponse depends crucially on the relative weight of the short-range component of 
disorder. Depending on the properties of disorder, the theory allows one to identify the temperature 
range within which the photoresponse is dominated by one of the mechanisms analyzed in the paper. 
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I. INTRODUCTION 

During the last decade many experimental and theoret- 
ical advances have been made in the field of nonequilib- 
rium magnetotransport in 2D electron systems. The re- 
search activity was triggered by the discovery of a number 
of beautiful nonequilibrium phenomena, the microwave- 
induced resistance oscillations (MIRO)-i 2 - in the first 
place, governed by the ratio lu/uj c of the radiation fre- 
quency and the cyclotron frequency cj g = |e| B/mc. Fur- 
ther experiments on the MIRO^"^ led to a spectacular 
observation of zero resistance states (ZRS) 3 -^^ which 
were later explained as a result of instability leading to 
the formation of current domains^ 4 - 

These discoveries stimulated an intense theoretical 
research 3 ^"— which has led to significant advances in un- 
derstanding the nonequilibrium transport in high Lan- 
dau levels (LLs). Initially, the MIRO were attributed to 
the "displacement' mechanism which accounts for spatial 
displacements of semiclassical electron orbits assosiated 
with the radiation-assisted scattering off disorderj 35 i 36 i 37 
The preferred direction of such displacements along the 
symmetry-breaking dc field oscillates as a function of 
uj/u> c due to a periodic modulation in the density of 
states (DOS) v{e) = v{e + fku c ) in a magnetic field. The 
displacement mechanism results in temperature indepen- 
dent MIRO with the phase and period as those observed 
in the experiments. 

The dominant contribution to the MIRO was later 
shown to be due to the "inelastic" mechanism associated 
with a radiation-induced changes in the energy distri- 
bution of electrons . 38 ! 39 The nonequilibrium part of the 
distribution oscillates with energy due to the oscillatory 
DOS and its amplitude is controlled by the inelastic re- 
laxation time ri n oc T -2 . 3 — As a consequence, the inelas- 
tic contribution to the MIRO provided an explanation 
of the observed strong temperature dependence of the 
MIRO. 

Further, a recent systematic study of the photore- 



sponse revealed two additional mechanisms of the MIRO, 
" quadrupole" and "photovoltaic." 4 ® In the quadrupole 
mechanism, the microwave radiation leads to excitation 
of the second angular harmonic of the distribution func- 
tion. The dc response in the resulting nonequilibrium 
state yields an oscillatory contribution to the Hall part 
of the photoconductivity tensor which violates Onsager 
symmetry^ In the photovoltaic mechanism, a combined 
action of the microwave and dc fields produces non-zero 
temporal harmonics of the distribution function. The re- 
sponse of this state to an ac electric field contributes to 
both the longitudinal and Hall dc photoconductivity. 

The theory developed in Refs. HHHO assumed that 
the only relevant source of disorder potential is that cre- 
ated by remote ionized donors which are separated from 
the plane of 2D electrons by a wide spacer of width 
strongly exceeding the Fermi wavelength, d 3> kp 1 . 
Such smooth disorder is characterized by scattering on 
small angles ~ l/{kpd), resulting in the long trans- 
port scattering time, r tr ~ Tqikpd) 2 ^> r q , and the 
backscattering time, tv oc r q exp(2fcj?d); here r q is the 
total disorder-induced quantum scattering time. For 
a smooth disorder potential, the inelastic contribution 
dominates over the displacement contribution for r- m > 
T q , the condition which is fulfilled in the whole temper- 
ature range where the MIRO were observed in the early 
experiments! 1 * 2 -^^ The quadrupole and photovoltaic 
mechanisms are the only ones yielding oscillatory correc- 
tions to the Hall part of the photoconductivity tensor, 
with the magnitude comparable to that of the displace- 
ment contribution to the diagonal part. 

In addition to the MIRO, 2D electrons exhibit nonequi- 
librium magnetooscillations when a strong current is 
applied to a sample in the absence of the microwave 
radiation i 65 i 66 ' 67 i 68 ' 69 These oscillations in the nonlin- 
ear resistivity, called hall-induced resistance oscillations 
(HIRO), were attributed to the resonant backscattering 
off disorder between LLs tilted by the electric fieldj 65 i 70 
Therefore, in contrast to the MIRO, the HIRO require 
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a sufficient amount of short-range disorder in order to 
be visible in experiment. This disorder potential is 
created by residual charged impurities in the vicinity 
of the 2D electron plane. Being inevitably present in 
quantum Hall structures, such impurities are believed to 
limit the mobility of ultra-high mobility samples. Re- 
cent experimental 2 - ! 29 ! 30 and theoretical 6 ^ study of the 
interplay of the MIRO and HIRO demonstrated that 
the displacement contribution to the MIRO is very sen- 
sitive to the nature of disorder. In fact, it is known 
that various transport phenomena in 2D electron sys- 
tems crucially depend on the type of disorder, includ- 
ing quantum magnetooscillations in the dc and ac trans- 
port, momentum-dependent conductivity, quasiclassical 
memory effects, commensurability oscillations in lateral 
superlattices, and interaction-induced quantum magne- 
toresistance (for a review, see Ref . l7ll) . Therefore, gen- 
eralization of the previous studies on the MIRO to the 
case of a generic disorder potential is highly desirable. 

Here we present a comprehensive analysis of the pho- 
toresponse, taking into account all four contributions to 
the MIRO» and using a generic-disorder model which is 
characterized by an arbitrary dependence of the elastic 
scattering rate 



n=— cx 



(1) 



on the scattering angle 8. We analyze the obtained re- 
sults for a realistic model of mixed disorder formed by 
superposition of a smooth random potential and short- 
range scatterers. The theory allows one to identify the 
temperature range where a specific mechanism dominates 
the photoresponse for different types of disorder. 

The paper is organized as follows. In Sec. [H] we present 
a kinetic equation approach to the problem. In Sec. IIIII 
we obtain the photoconductivity tensor which is analyzed 
in Sec. IIVI for different types of disorder. In Sec. [V] we 
disccuss possible temperature regimes of the MIRO tak- 
ing into account the interaction-induced broadening of 
Landau levels. Main findings are summarized in Sec. I VII 



II. FORMALISM 

To begin with, we briefly outline the main steps in 
setting up the formalism. We consider 2D electrons in a 
classically strong magnetic field B in the presence of a 
random potential, a dc electric field, 



Edc — (E x , E y ), 
and a microwave field 

E u (t) = ^^ S± Re[e±e^ t ] 



(2) 
(3) 



where \ 2e± = e x ± ie y , and s± parameterize polariza- 
tion of the microwaves (s^_ + s_ = 1). The main param- 



eters in the problem are related to each other as follows: 

E F > T, U , UJ C , T q _1 > Tt," 1 , 

where sf is the Fermi energy and we use ks = ft = 1. 
We adopt the approac h 37 ! 39 ' 40 ! 62 to the problem, based 
on the quantum Boltzmann equation (QBE) for the semi- 
classical distribution function f(e,tp,t) at higher LLs, 



(d t + w c S„)/ + rr^if) - f T ) = St{/} . 



(4) 



In the inelastic collision integral oc t-~ , /t(e) is the 
Fermi-Dirac distribution. The angular brackets (...) de- 
note averaging over the direction ip of the kinematic mo- 
mentum. 

The QBE (|3|) allows us to treat the interplay of the dis- 
order, the Landau quantization, and the external fields, 
which are all included into the impurity collision integral 
St{/}. The field-dependent collision integral appears as a 
consequence of transition to a moving coordinate frame, 
r — > r — £ t , where 

and i xy — —i yx — 1. In the new frame, the electric field 
is absent, but the impurities are moving and can change 
the energy of electrons. The effects of generic disorder 
and external fields are encoded in the integral operator 62 



dip' e ipp ( Tl *' - 'V)^ t 2 - ^ t i) 



2tt 



Tip-pi 



(6) 

where pf — vtivf is the Fermi momentum, the unit vector 
n v = (cos ip, sin ip), and F(tp) is an arbitrary function. In 

terms of /C, the scattering integral in the time represen- 
tation is given by the expression 



st{/} 



dU 



^t 2 ..t 1 {g^ 2 -t 3 ft a ,t 1 - /t2,ts5t 3 -ti} 



(J) 



We are going to explore the regime of overlapping LLs, 
^cT q <C 1, when the spectral functions 

= K-J - Kt2 - ti)/2 - XS(t 2 - tl - t B ), (8) 

depend on a single parameter 

A = exp(-7r/w c r q ) <C 1 (9) 

and are insensitive to the external electric fields. The 
corresponding density of states (DOS) is 



v{e) = vq[1 — 2Acos(ets)], 



(10) 



where ts — 2tt/|o; c | is the cyclotron period and vq = 
m/2-K is the DOS at B = per spin degree of freedom. 
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The Wigner transform (Fourier transform with respect 
to t_ = t 2 — ti) of the scattering integral, Eq. ([7]), defines 
the QBE for the distribution function f(e,ip,t), which 
can be approximated by 



/ = h + (d s f T ){4ti - 2ARe[0 1 ex V (iet B )} + 2A 2 ^ 2) } 



where 



4>q(<p) = ev F E dc ■ n v /uj c Z T tr 



(11) 



(12) 



is the classical part leading to the Drude expression for 
the current §TE\i , and t = (t 2 + t\)/2. By the symmetry 
of the kernel ©, 



(13) 



includes even angular harmonics of the distribution func- 
tion and is governed by the real part of the kernel, Eq. ([6]): 



(d t + u c d v )<f>± + ' 



fa 



di 



jRe/C t2 , tl {l} 



(14) 



t 2 =t 1 +t I 



The dc electric current, 



2ev F / d£v(e)(n v f(£,(p,t)) -2eu e F d t C(t) , (15) 



is determined by the first angular harmonic of the distri- 
bution function, Eq. (jlip , which is present in Re^i , <pQ D ^ > 

(2) 

and </>q . The bar denotes time averaging over the period 
of the microwave field in the steady state. 

Our aim is to calculate the current ([15)) for generic 
disorder |T]) to order E 2 Ed c for temperatures T » ui,ui c . 
This regime is the most relevant to experiments on the 
MIRO and allows for a reliable comparison between the 
theory and experiment. In the absence of quantum cor- 
rections (A = 0) the current is given by the Drude for- 
mula, 



j D = 2a D {1 - u) c T tI £)E d c, 



(16) 



where <td = e 2 vov F /2ui c 2 Tt I . Neither strong dc field 
nor microwaves modify this result as long as A = 
and the energy dependence of elastic scattering (weak 
at uj c <C £p) is neglected4£ All quantum corrections to 
v{e)f(e) of first order in A oscillate with energy and 
are exponentially suppressed at 2tt 2 T jiO c 3> 1 similar 
to Shubnikov-de Haas oscillations. The leading quan- 
tum corrections which survive the temperature smearing 
at T > a) c are of order A 2 and can be written using 
Eqs. flTTJ), (T5J) as 



j =i D - 4\ 2 ev F v Q (ReJ,lmJ) T 



(17) 



with 



J = 



(>"r I - - 

— d t _lC t 2jtl {l} - 0j_/C t2jtl {l} 



lCt 2 ,t 1 {(^±}) 



(18) 



III. PHOTOCONDUCTIVITY TENSOR 

We show that the current J [Eq. (fT8|) ] depends in an 
essential way on the properties of disorder; specifically, 
on whether it scatters isotropically (short-range disorder) 
or primarily on small angles (long-range disorder). A 
straightforward calculation of J to the first order in Ed c 
and to the second order in E w using the solution <f>±_ of 
Eq. ([TJ} gives the following expression for the current 
(HZ}: 



j =j D +4,a D X 2 E dc + a ( P h '>E dc 
h s — h a d s d a 



(19) 
(20) 



The terms of different symmetry d Sta and h s _ a in the pho- 
toconductivity tensor cr(P h ) are given below for arbitrary 
strength of the short-range and smooth components of 
disorder. The diagonal part of the photoconductivity 
tensor consists of the isotropic, d s , and anisotropic, d a , 
components, where 

d s = d{ A *> + df> + d{ D \ (21a) 

d{ A > 



A 2 [sin 2 w + w sin 2w] £ 2 



4 s ' = A>sin2 W ]^^£ : 

7"tr 



2 



di 1 



-A 2 [w sin 2 w] 



2 { 2r tr f~ 2 

LUT tv LO ± 2UJ C 



and 



d a = d^+di D \ 



(21b) 



(21c) 



(21d) 



(21e) 



d{ A) = -A 2 [sin 2 w + wsin2w]^-£ + £_ , (21f) 

2r* 



AD) = A 2r„..„ : „2 



w sm~ w\ £+£- 

LOTtr 



(21g) 



In the off-diagonal part, h a represents the microwave- 
induced correction to the Hall conductivity, while h s is 
the anisotropic contribution to the dissipative conductiv- 
ity (violating the Onsager symmetry) 40 ' 64 

h a = h ( f ) = -\ 2 [wsm2w}(T tI /oj c f 2 )£ + £- , (22a) 
K = h a D) = -A 2 [2 sin 2 w + w sin 2w] 



2>4(i 



TtrT 2 



OJTti CO ± 2iO t 



(22b) 



All microwave-induced corrections (|2"Tj) and (|22"|) show 
oscillations with the ratio 

LOtB 7T07 

W ~ ~2~ ~ IwJ 



(23) 



and are of the second order in the microwave field E c 
measured in dimensionless units 

ev F E w 



£ 



± 



oj(u> ± uj c ) 



(24) 
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In Eqs. (|21j) and (p2|) . the superscripts denote the 
contributions from the displacement (A), inelastic (B), 
quadrupole (C) and photovoltaic (D) mechanisms.™ The 
displacement contribution, (A), comes from the first term 
in Eq. (fT8|) . In the absense of the microwave radia- 
tion, E u = 0, this term produces the quantum correction 
to the Drude conductivity given by the second term in 
Eq. ([15]). At order E Ac E^, the first term in Eq. (Tn§ ac- 
counts for the microwave-assisted disorder scattering in 
the presence of the symmetry-breaking dc field. There- 
fore, in the displacement mechanism the radiation di- 
rectly affects the first angular harmonic of the distribu- 
tion function, and thus the dc current, see Eqs. (|21b[) 
and (|21fj) . Mechanisms (B), (C), and (D) are related to 
the microwave excitation of even angular harmonics of 
the distribution function governed by Eq. (1141) . Their 
contribution to the current appears at the same order 
Ed c E^ and is described by two last terms in Eq. (fl~8|) . 
The microwave-induced changes in the time-independent 
zero angular harmonics of the distribution function are 
controlled by the inelastic relaxation, and lead to the in- 
elastic contribution, (B), Eq. (|21c|) . In the quadrupole 
mechanism, (C), the microwave radiation excites the 
second angular harmonic of the distribution function, 
oc E 2 cos 2(f>, which also contributes to the dc cur- 
rent at order Ea c E 2 after substitution into Eq. (fT8j) with 
JC oc Edc, see Eq. (|22aj) . In the photovoltaic mechanism, 
(D), a combined action of the microwave and dc fields in 
Eq. (|14j) produces both zero and second angular harmon- 
ics of the distribution function which oscillate in time 
with the microwave frequency. The ac response in the 
resulting state [Eq. (fig)) with JC oc E u and <j)±_ oc E^E^] 
contributes to the longitudinal and Hall parts of cn ph ), 
see Eqs. (|2Td) , fl2Tgl ), and (p2b|) . 

The photoconductivity tensor ([20]) depends on the in- 
elastic scattering rate -q" 1 which enters the inelastic con- 
tribution (|21c| and on four different disorder scattering 
rates. The latter are expressed in terms of the angular 
harmonics r„ of the disorder scattering rate, Eq. |T]), as 



(25a) 
(25b) 
(25c) 
(25d) 



We emphasize that as compared to the case of smooth 
disorder, where the rates t _1 and t^ 1 fully parametrize 
the photoconductivity,— two additional rates r" 1 and 
f _1 are required in general situation. The rate t~ con- 
trols the magnitude of the displacement contributions 
(|21b|) and (|2Tf]l . reproducing the result of Ref. El. Note 
that, in view of <j]j oc t^ 1 , the product ajjdi A ^ oc t^ 1 
entering Eq. (|20p is actually independent of r tr . The rate 
f _1 enters the quadrupole and photovoltaic contributions 



1 


1 






Tq 


! 

TO 






1 


1 


1 




Ttr 


TO 


7 

Tl 




1 


3 


4 


1 


t* 


TO 


Tl 


T 2 ' 


1 


1 


1 


1 


f 


Tl 


~ 2^ 


"2^ 



(|21d [) , plg] )-(|22b |) . The relation between the scattering 
rates (|25p strongly depends on the type of disorder po- 
tential in the plane of 2D electrons (short-range vs. long- 
range). We discuss this dependence in the following sec- 
tion for a realistic model of the disorder. 



IV. DEPENDENCE OF MIRO ON TYPE OF 
DISORDER 

In this section we study the relative magnitude of the con- 
tributions of mechanisms (A)-(D) for different types of 
disorder. Using Eqs. (f21~j) . (|22|) and ([25]) we demonstrate 
that the MIRO amplitude is very sensitive to the details 
of the disorder potential and may provide a valuable in- 
formation about the disorder which can not be extracted 
from the conventional transport measurements of the mo- 
bility and the Shubnikov-de Haas (SdH) oscillations. 74 



A. Mixed-disorder model 

The major source of elastic scattering in high-mobility 
Hall structures is a smooth random potential created by 
remote donors that are separated by a large spacer of 
width d ^s> kp 1 from the 2DEG plane. The disorder 
model which takes into account only such remote charged 
impurities was used in a number of theoretical works on 
the MIR 37 ' 39 i 4 ° since it allows one to account for the 
experimentally relevant small-angle scattering condition 



T q <C Tt r and to consider in an unambiguous way the re- 
gion of magnetic fields r^ 1 <C w c < T q 1 ■ Here we imple- 
ment a more realistic 7 ^ "two-component" (or "mixed") 
model of disorder— which includes strong short-range 
scatterers in addition to the smooth potential. 

In terms of the angular harmonics of the elastic scat- 
tering rate ([1]), the mixed disorder model is formulated 
as 



O 62,70 



1 



3n,0 
T s h 



1 



1 



1 + xn 2 



(26) 



Here the isotropic scattering rate tJZ characterizes the 
short-range disorder created by residual impurities lo- 
cated at or near the interface, while the smooth part 
is defined by the total (quantum) rate t - ^ 1 and by the 
parameter \ — (kpd)~ 2 <C 1 giving a characteristic scat- 
tering angle 9 ~ <§; 1. 

Since in high- mobility samples r tr 3> r q , we imply that 
tT <C ^~sm an d the smooth component dominates in the 
quantum relaxation rate 

11 1 1 

Tq T T sh T sm 



(27) 



At the same time, the relative weight of the sharp and 
smooth components in the transport relaxation rate 



1 



Ttr 



1 



Tsh 



1 



X 



1 + X 



(28) 
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can be arbitrary. 



B. Photoconductivity for smooth disorder 

It is useful to recall the results of Ref. HO for the case of 
smooth disorder and to see how they are obtained from 
Eqs. (EH), (HH) and (JUJ. Putting t^ 1 = (no short- 
range disorder), the scattering rates that enter Eq. 
can be rewritten as 



1 


1 


(29a) 


Tq 




1 


x_ 

5 


(29b) 


7tr 






1 


12x 2 


(29c) 




3 






1 

f 


1 
T tr ' 


(29d) 



As a result, in the case of smooth disorder the ratio 
of the inelastic contribution (|21c|) to the displacement 
contribution (|21bp is 



4Ti n T* 



3t q 



'sh 



0. 



(30) 



This relation follows from Eqs. (|21bj) and (I21cj) if one dis- 
regards the term oc sin 2 w which is small at w = ituj /ui c ^> 
1. 

For typical conditions of the MIRO experiments, the 
parameters entering equations Eqs. (f2"Tj) . (|2"2"j) and (j2"9"|) 
can be estimated as w c ~ w ~ r q~ 1 an d X = T q /r tr ~ 
1CP 2 . Using the estimate T m ~ £f/T 2 for the inelas- 
tic scattering time, we obtain Tj n ~ r tr at T = IK. 
Therefore, the ratio (f3"0"|) remains small up to T ~ 10-ftT 
(where the MIRO are strongly suppressed). The mag- 
nitude of the oscillations produced by the quadrupole 
and photovolatic mechanisms is of the order of the dis- 
placement contribution, as they contain a small prefactor 

Kr)" 1 ~ (cJ c T tr ) _1 ~ Tq/r tr . 

We come to the conclusion that, in the case of smooth 
disorder, the photoconductivity tensor is dominated by 
inelastic mechanism (B). It produces the isotropic di- 
agonal contribution (|21cp proportional to r; n ~ ep/T 2 . 
Other mechanisms produce weak anisotropy of the photo- 
conductivity tensor [Eqs. (|21fj) . ( |21g[ ), (|22a[) ] and govern 
the Hall part of the photoconductivity, (|22b[) . 



C. Photoconductivity for mixed disorder 

Now we return to the mixed-disorder model (|26[) with 
nonzero r^ 1 and analyze the dependence of various con- 
tributions to the photoconductivity tensor, Eq. (|20p . on 
the weight x — r tr /T s h of the sharp component of disor- 
der in the transport relaxation rate (|28p . In this analysis, 
we fix both r tr , Eq. l|2"8|) , and r q , Eq. I|27p. and use as a 



parameter the ratio a = r tr /T q which can be extracted 
from the experiment^ 

Inspection of Eqs. (f2"Tj) . (|22p shows that only the in- 
elastic contribution (I21cp is independent of the type of 
disorder at fixed Tt r and r q , 



d{ BS> cc ri n /Tt r = const(x) . 



The displacement contributions d 



(A) 



(31) 
and da , 



Eqs. (|21bp and (|21fp . are proportinal to the ratio 

ax — bx + 4 



111 

T* 



4x + 3 ' 



(32) 



which is illustrated in Fig. QJi for several values of the pa- 
rameter a — r tr /r q . Since in high-mobility structures the 
parameter a is always large, Eq. (|32p demonstrates that 
the displacement contribution is parametrically enhanced 
in the case of short-range disorder, \ x =i/t*\ x =o = 4ct > 
1. 




0.2 0.4 0.6 

ttAsh 



FIG. 1: Dependence of the relaxation rates r, _1 (a) and 
Tu/f 2 (b) on the type of disorder (parametrized by the weight 
Ttr/Vsh of the short-range component of disorder in the trans- 
port relaxation rate) for fixed r tr and r q for three values of 
a = 100,30, 10. 

For the ratio (rt r /f) 2 which enters the photovoltaic 
and quadrupole contributions (|21dp ~ (|22bp . we obtain the 
dependence illustrated in Fig. []J>, 



2a — 3ax + 7x — 6 
2a - 8x + 6 



(33) 



Unlike the displacement contribution, Fig. QJi, the pho- 
tovoltaic and quadrupole contributions are at their max- 
imum in the absence of the sharp component, x = 0, 
and do not change parametrically in the opposite limit: 
f~ 2 \ x =i/f~ 2 \ x =o = 1/4 for a 3> 1. However, the two 
contributions change significantly as a function of x be- 
tween x = and x = 1 since at x = (2a — 6) /(3a — 7) 
the effective rate f _1 changes sign. 

Summarizing, for a sufficiently large concentration of 
short-range scatterers, i. e. x ~ 1 and r sri ~ r tr , the 
ratio Ttr/r» controlling the displacement contribution, 
Eqs. (|21bp and (f21~fj) is of order unity (see Eq. (JS2) and 
Fig. For typical experimental parameters, the ratio 
Tm/7"tr entering the inelastic contribution (|21c|) is also of 
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order unity at T ~ IK (see Sec. IIVBI) . Therefore, un- 
like the case of smooth disorder discussed in Sec. IIVBI in 
the present scenario of strong short-range disorder both 
mechanisms produce comparable contributions to the di- 
agonal isotropic component d s of the photoconductivity 
(f20|) . At the same time, the photovoltaic and quadrupole 

contributions di D J , h s and h a still contain the small fac- 
tor (wr tr ) _1 ~ (w c r tr ) _1 which is of order r q /r tr for 

w c r q ~ 1. Consequently, the contributions di^J to the 
diagonal part of the photoconductivity tensor can be ne- 
glected. The terms h s and h a are also small; however, 
they remain the only microwave-induced contributions 
to the off-diagonal part of the photoconductivity tensor 
and, therefore, are important. 

Since the displacement and inelastic mechanisms are 
equally important in the diagonal photoresponse in the 
case of strong short-range disorder, x *~ 1, it is natural 
to discuss possible ways to separate their contributions 
experimentally. As we argue below, the difference in the 
temperature dependence of the two contributions can be 
used to extract the contributions separately. In Sec. [V] 
the temperature dependence of various contributions is 
analyzed in detail. 



TEMPERATURE DEPENDENCE OF MIRO 
FOR SMOOTH AND MIXED DISORDER 



Temperature dependence of MIRO for 
disorder-broadened Landau levels 



In the photoconductivity tensor, as given by Eqs. (|2U|) . 
(|2"Tj) and (|2"2"|) . the only T-dependent quantity is the in- 
elastic scattering time Ti n entering the inelastic contri- 
bution (|21c[) . According to the analysis in Ref. [H, the 
inelastic relaxation of the e-oscillations in the distribu- 
tion function is dominated by electron-electron (e-e) col- 
lisions. Under conditions of the MIRO experiments, the 
e-e collision rate for an electron at energy e counted from 
the Fermi level Ef is 



1 



e 2 + tt 2 T 2 



r ee (e,T) inei 



In 



kvf 



max{T, cj c (cj c T tT ) 1 / 2 }' 



(34) 



with the inverse screening length k = 4ire 2 Vo.— The ef- 
fective inelastic scattering time entering (|21c[) is given by 
the thermal average of the out-scattering time, 

^ = {t 00 (£, T)} T ^ 0.822 r cc (0, T), (35) 

where the notation 



{Feh 



deF £ d £ f T (e) 



(36) 



has been introduced. 

The above results are applicable as long as the disor- 
der alone controls both the momentum relaxation and 



the energy dependence of the DOS. Namely, when the 
contributions of the inelastic scattering mechanisms to 
the momentum relaxation rate, T tr _1 , and the quantum 
scattering rate, Tq" 1 , are negligible. This is usually the 
case at T ~ IK, which is a typical temperature in MIRO 
experiments. 



B. Effect of interactions on the density of states 

At higher temperatures, the momentum relaxation 
is modified by the phonon effects, while the DOS 
is influenced by both the e-e and elect ron-phonon 
interaction . 27 i 32 ' 42 i 50 In the limit of overlapping LLs, the 
effect of the interactions on the DOS amounts to substi- 
tuting A in Eqs. (fl"5|) and (fTTj) with 



A(e,T) = exp - 



^ 1 +^c 1 (£,r) + r c l 1 ph ( £ ,T) 



(37) 

The phonon effects will be considered elsewhere^ 2 - Here 
we discuss the modification of the temperature depen- 
dence of the photoconductivity tensor (|20|) . (|2"Tj) . (f22|) 
caused by the effect of the e-e interaction on the DOS 
(in recent experiment s 27 ' 32 the e-e interaction dominates 
the effect on the DOS). As shown below, the magnitude 
of the above effect is controlled by the parameter 



2tt 



(38) 



C. Interaction-induced exponential T— decay of the 
MIRO 



In order to account for the interaction-induced variation 
(f3"T)) of the amplitude of the DOS oscillations (which is 
slow on the scale lo c <C T), we make the replacement 



\ 2 = {\ 2 } T ^f 1 (T)={~\ 2 ( S) T)} 1 



(39) 



everywhere in Eqs. (121]) and ([22]) except for the inelas- 
Eq. 

replacement is, see Eq. (|55l) . 



tic contribution di B \ Eq. (I21cl) . In the latter case, the 



A 2 — 

Ttr 



/ 2 (T)= A 2 (e,T) 



r ee (s,T) 



Tu- 



rn 



At 7 -C 1, the interaction- induced modification of the 
DOS is small, yielding 



h = A 2 (l- 7 /3) 
f2 = \ 2 (- 



7< 1. 



2tt 



Ttr 0J c T tl 



7 < 1 



(41) 
(42) 



The correction to the inelastic contribution d s is tem- 
perature independent, while all other contributions ac- 
quire a weak temperature dependence through 7 cx T 2 . 



7 



In the opposite limit 7> 1, 



ft 
h 



A 2 
A 2 



4 7 V2 ' 

^ 5 / 2 e^ 
w c r tr 



2 7 3 / 2 



7> 1. 
7 > 1. 



(43) 
(44) 



In this high temperature limit, all the quantum effects 
oc A 2 that survive the thermal averaging at 2tt 2 T/ui c ^> 1 
are exponentially suppressed by the inelastic e-e scatter- 
ing. 

Interestingly, the SdH oscillations behave just in the 
opposite way: They are exponentially suppressed by 
thermal smearing at 2tt 2 T/uj c > 1 but they are not 
influence d 77 ' 78 ' 79 by the interactions that lead to the 
dependence of t^} on e in the form t^} cx e 2 + 
7r 2 T 2 . Indeed, for such a specific form of r" 1 entering 




(b) 



3 4 
T, K 



3 

T, K 



Eq. (|37p , the thermal average produces the unperturbed 

result, {A(e, T) cos{2-ke/uj c )}t — * A{cos(27T£/w c )}t for 
T lu c , because in this case the energy-dependent 

term in r" 1 effectively shrinks the range of energy in- to ^ mo ^, h dis ™ der ^i in r stance ' , T * r /T s , h = , and Ttr / r ' 
tegration and increases {cos(27re/ti; c )}T to compensate 
exp(-7r/w c r co (0,T)). 



FIG. 2: Solid lines: T-dependence of the magnitude of 
the inelastic contribution ipic|l . 2{r co (e, T)A 2 (e, T)} T /r tr A 2 , 
for 7iii|t=ik — 2r tr . Dashed lines: T-dependence of 
the magnitude of the displacement contribution ()21b|l . 
Tt r {A 2 (e, T)}t/2t» A 2 , for r tr /2r» = 1 (which corresponds to 
Ttf/Tsh = 0.6 — 0.7 in Fig. [TJ. Dotted lines: T-dependence 
of the magnitude of the displacement contribution (I21b[l . 
Ttr{A 2 (e,T)}T/2r»A 2 , for r tr /2r» = 0.1 [which corresponds 

50 

in Fig. [T] see Eq. (|32p]. Left and right panels correspond to 
(a) weak effect of interactions on the DOS at T = 1 K, 71 = 
2n/uj c r cc (eF, 1 K) = 0.01, and to (b) strong effect, 71 = 0.1. 



D. Temperature regimes of the MIRO for smooth 
and sharp disorder 

It follows from the above arguments that different tem- 
perature dependences of magnetoresistance are expected 
for samples with smooth and sharp disorder. In the for- 
mer case Ttr/V* <§C 1, and the inelastic contribution (|21cp 
is bigger than all other contributions by a factor of Ti n /T q , 
see Eq. (l30l) . This ratio is estimated as Ti n /r q ~ 10 — 10 2 
at T = 1 K, making the inelastic mechanism dominant 
in the whole range of temperatures where the MIRO are 
observed. The temperature dependence of the diagonal 
isotropic component d s of the photoconductivity tensor 
(|2"Tj)) in the case of smooth disorder is then given by the 
function f 2 , Eqs. 1(10]). (|4"2"j) . and (|3I|). which is shown in 
Fig. [2] by the solid line. Specifically, the T~ 2 behavior 
(|42| at low T crosses over to the exponential decay (|44"]) 
at high T. 

The situation changes in the case of strong short-range 
component of disorder, T tr /r s h ~ 1. Since T tr /r* ~ 1 in 
this case, the displacement contribution (|21b[) to the di- 
agonal isotropic component d s of the photoconductivity 
tensor ([2D]) starts to dominate over the inelastic contri- 
bution (|21c[) at a certain characteristic temperature T\ 
at which the MIRO are still strong. The temperature T\ 
in Kelvin is given by the relation 

Ti = lKx — \J t» Tj n (l K), (45) 

Ttr 

see Eqs. ([30|l . ([35]). and Experimentally, T x is of 

order 1 K for the case of strong sharp component of dis- 
order since both T»/T tr and Tj n (l K)/rt r are of order unity. 



Another characteristic temperature T2 is given by the re- 
lation 7(T 2 ) = 1, see Eqs. (J35J and J3SJ, 

T2 = y/a^K) - 1 K x 0.44V^ c r in (l K). (46) 

It separates the temperature region T > Ti where the 
LLs broadening induced by the e-e interaction becomes 
essential. 

In the case T\ -C T2 the diagonal symmetric part 
d s of the phoconductivity (|2U)l has different tempera- 
ture dependence in three temperature intervals sepa- 
rated by Tx and T 2 . For T < Tj the inelastic mecha- 
nism is dominant and, as a result, d s ~ di B%> cx T -2 . 
In the interval T\ < T < T 2 , the inelastic contribu- 
tion becomes smaller than the displacement contribution 
leading to a temperature independent magnetoresistance 
d s ~ ~ const(T). Finally, for T > T 2 all the con- 
tributions are suppressed exponentially, see Eq. (|4ip and 

Fig.m 

Clearly, the T-independent interval doesn't exist when 
T2 is of order or smaller than T±. This situation was 
probably realized in the recent experiment^. The pho- 
toresponse is then expected to cross over from the inverse 
quadratic dependence d s oc T~ 2 for T < T 2 to a much 
faster exponential decay for T > T 2 . In this case, it is 
difficult to separate the displacement and inelastic contri- 
butions using the temperature dependence of the MIRO 
only, since in the region where the two contributions have 
comparable magnitudes, their temperature dependence 
is almost identical (see Fig. Wp)- Therefore, for Tx > T 2 
some additional measurements aimed at determination 
of, for instance, the anisotropic diagonal component d a 
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of the photoconductivity (|20[) are desirable. Since d a , 
Eq. (|21fj) . is produced solely by the displacement mech- 
anism (the photovoltaic diagonal contributions and 

d{ D ^ are negligible for Tt r /r s h ~ 1), such a measurement 
would be a reliable tool to estimate the significance of the 
displacement mechanism at a given temperature. Alter- 
natively, for a fixed polarization of the microwaves, one 
can analyze the photoresponse in the vicinity of "odd 
nodes" lo/lo c = n + 1/2, where the displacement con- 
tribution (|21b|) is non-zero, sin (ttu>/uj c ) = 1, while the 
inelastic contribution oc sin(27ru;/a; c ) = vanishes. 

VI. CONCLUSION 

In conclusion, we studied the microwave-induced magne- 
tooscillations (MIRO) in the conductivity of a 2D elec- 
tron gas with mixed disorder containing both a long- 
range random potential and short-range scatterers. We 
calculated all contributions to the current that are lin- 
ear in the applied dc electric field and quadratic in the 
microwave field. It is found that the relation between dif- 
ferent contributions strongly depends on the strength of 
wide-angle scattering off disorder and the sample temper- 
ature. In general, the microwave radiation modifies both 
the diagonal and off-diagonal components of the conduc- 
tivity tensor, including the Hall component. Study of 
anisotropy of the current may identify each contribution 
separately. The dominant contribution at low temper- 
atures is the inelastic contribution, which is always di- 
agonal and isotropic. At higher temperatures and in 
the presence of wide-angle scattering, the displacement 
mechanism may become larger than the inelastic contri- 
bution. At high temperatures the quantum corrections 
to the conductivity are exponentially suppressed by the 
electron-electron interaction. 

We believe that the possibility of measuring two types 
of nonequilibrium magnetooscillations — the MIRO and 
the HIRO — in the same sample has opened new avenues 
for the characterization of 2D electron systems. In par- 
ticular, measuring both of them could unambiguously re- 
veal the two-scale structure of disorder characteristic of 



high-mobility samples. Indeed, the relative weights of 
the smooth component of disorder (produced by remote 
impurities separated from the electron gas by a wide 
spacer) and the short-range component (produced by 
residual impurities sitting near or at the interface) mani- 
fest themselves differently in the various relaxation times 
that could be extracted from the experimental data. As 
we have demonstrated in this paper, the MIRO are sen- 
sitive to the correlation properties of disorder. Specif- 
ically, their amplitude and shape are parametrized by 
the inelastic relaxation time and the disorder-induced 
relaxation times for the zeroth, first, and second angu- 
lar harmonics of the electron distribution function, with 
the displacement contribution to the MIRO — strongly 
enhanced in the presence of short-range disorder — being 
parametrized by the time r*. Therefore, in addition to 
the routinely measured TsdHCr^ and Tt r , one can: (i) 
extract the time r q from the exponential dependence of 
the amplitude of the oscillations on B, both the MIRO 
and the HIRO; (ii) extract the time t* from the T and 
u> dependences of the MIRO, as well as from the ampli- 
tude of the HIRO. As for the inelastic time, it can be 
extracted in a variety of ways: from the power-law T de- 
pendence of the MIRO amplitude at lower T, from the 
MIRO and/or HIRO exponential behavior as a function 
of T at higher T, as well as from the behavior of the non- 
linear dc resistance as a function of the applied dc electric 
fields Importantly, the inelastic and displacement con- 
tributions to the MIRO can then be unambiguously sepa- 
rated from each other. One interesting possibility worthy 
of experimental investigation is to extract all the relevant 
parameters in the absence of microwaves and use them 
to describe the experimental data on the MIRO without 
any free parameters. 
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